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1.0 SUMMARY 


This is the final report for the Control Technology Aeroassisted 
Orbital Transfer Vehicle study contract (NAS8-37358) . 

The long range goal of the Aeroassisted Orbital Transfer Vehicle 
(AO TV) program is to furnish NASA with a reusable orbital transfer 
vehicle with a larger payload capability, achieved by replacing a 
propulsive orbit transfer maneuver with an aeroassisted maneuver. 
The focus of this contract has been to develop the control 
technology required to identify the sophistication required for 
the AO TV control system. Our primary control objective has been 
to develop control techniques which minimize the exit condition 
errors while minimizing the control effector energy usage. The 
control technology also promotes low cost, high reliability, 
serviceability, refueling, recalibration, autonomy, and 
space-basing. 

The primary design challenge for the AOTV mission is to minimize 
the fuel used for atmospheric maneuvering and orbit establishment 
following the aerobraking maneuver. (Fuel used for orbit 
establishment is estimated by assuming impulsive velocity changes, 
referred to as AV burns.) This can be accomplished with a 
combination of effective guidance and control algorithms. The 
HYPAS algorithm used in this study made it necessary to use large 
thrusters to give the vehicle a response fast enough to follow the 
quickly changing large angle guidance commands. 
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An angle-of-attack and bank angle command control technique has 
been developed consistent with the AOTV guidance approach. Both 
on-off thrusters and propo.rtional thruster systems have been 
analyzed for the control effectors. For the proportional 
thrusters, a classical linear control design was synthesized. For 
the conventional on-off thrusters, a phase-plane control algorithm 
was developed. Both designs furnish adequate control with the 
proportional thruster technique using less energy. An 
angle-of-attack adaptive controller was used to minimize the RCS 
fuel usage due to payload center of gravity uncertainties. A 
semi-adaptive gain computer was developed to minimize the mission 
data load. 

The control techniques were verified using a detailed 6D0F 
simulation. The effect of parameter uncertainty on mission 
performance and guidance/control performance was examined. 
Uncertainties in the entry state, mass properties, atmosphere, 
aerodynamics, and sensor data were evaluated. 

The vehicle performed well for most parameter variations, being 
the most sensitive to entry state errors and least sensitive to 
sensor and aerodynamic errors, suggesting that the system is very 
sensitive to the guidance and not as sensitive to uncertainties in 
the dynamics. This turned out to be the case as the primary 
difficulty in controlling the AOTV in this study was getting 
enough control authority to effectively follow the HYPAS guidance 
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commands . 


Several advanced control concepts are presented in the report. 
Advanced technology for the AOTV includes semi-adaptive and fully 
adaptive control techniques, as well as control designs which 
minimize cost and increase reliability. 

A nonlinear six-degree-of-freedom (6D0F) simulation was used to 
evaluate the attitude control techniques. The simulation included 
the Hybrid Predictive Aerobraking Scheme (HYPAS) guidance 
algorithm. Global Reference Atmosphere Model (GRAM) atmosphere 
models, and full nonlinear equations of motion including orbital 
effects .- 

A three-degree-o f- f reedom (3D0F) simulation was used for 
comparison purposes. The 3D0F simulation included the same 
guidance algorithm but only point mass dynamics and no control 
system. Outputs of the two simulations were compared to validate 
correct implementation of the guidance algorithm and equations of 
motion. 
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2.0 CONTROL CONCEPT 


The control concept for the AOTV was developed based on the 
mission requirements, the dynamic properties of the shaped-brake 
vehicle, and the operation of the guidance algorithm^ The 
following sections describe each of these elements and the 
resulting attitude control system. 

2 . 1 Mission 

Several AOTV concepts, including a ballute-type vehicle, a 
lifting-brake vehicle, and a shaped-brake vehicle, have been 
considered in recent studies, see Figure 2.1-1. The shaped-brake 
vehicle was used in this study. 

The vehicles are used to transfer payloads between geosynchronous 
(GEO) and low earth orbit (LEO) . A pass through the Earth's 
atmosphere, referred to as the "aeropass," is used to "brake" the 
vehicle for insertion into the low earth orbit. The aeropass 
reduces the amount of fuel required to slow the vehicle. Figure 
2.1-2 shows a transfer from GEO to LEO. Following the aeropass, 
three propulsive burns are used to establish the correct low earth 
orbit. These burns are also shown on Figure 2.1-2. A typical 
altitude versus velocity profile is shown in Figure 2.1-3. Notice 
that the vehicle exits the atmosphere (roughly 400,000 feet) 9500 
feet per second slower than the entrance velocity. Figures 2.1-4 
and 2.1-5 show the dynamic pressure and the bank angle versus. time 
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angle versus time profiles for the nominal mission, respectively. 


2 . 2 AO TV Shaped-Brake Characteristics 

The shaped-brake concept has the property that the lift vector can 
be directed by rolling the vehicle around the velocity vector. 
Thus, the vehicle can be used to pull into of out of the orbital 
plane as well as up or down, see Figure 2.2-1. The shaped-brake 
vehicle typically flies at large angles of attack (70° - 80°), 
meaning that the rolling motion about the velocity vector is 
principally a yawing motion in the body axes. The shaped-brake 
vehicle and the aerodynamic and control coordinate systems are 
shown in Figure 2.2-1. The aerodynamic characteristics of the 
shaped-brake are shown in Figure 2.2-2. For the vehicle used in 
this study, the L/D ratio was 0.2. 

2.3 AOTV Performance Corridor 

The shaped-brake AOTV has a lifting characteristic that can be 
used to control the vehicle trajectory and obtain a desirable exit 
velocity and inclination. The L/D ratio defines a region of 
successful mission operation, referred to as the control corridor, 
as illustrated in Figure 2.3-1. The name control corridor comes 
from the fact that within this region, the AOTV can be 
successfully controlled and guided to the exit window by directing 
the lift vector. Bigger L/D ratios correspond to wider corridors. 
If the trajectory carries the vehicle outside of the corridor, the 
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AO TV will not achieve the desired exit window. The upper corridor 
boundary is characterized by a continuous downward pointing lift 
vector, which is used to pull the vehicle into the atmosphere in 
order to decelerate to an acceptable exit speed. Conversely, the 
lower corridor boundary is characterized by a continuous upward 
pointing lift vector, used to avoid excessive speed depletion. 
Figure 2.3-2 shows the relationship between the total orbit 
insertion AV burn and the control corridor. Inside the corridor, a 
minimum AV burn can be achieved rotating the lift vector such that 
the desired trajectory is obtained. 

The control corridor width is closely related to the lift 
characteristic of the vehicle; with a high L/D ratio, the AOTV has 
more control authority and is more capable of pulling in and out 
of the atmosphere. Conversely, at a low L/D ratio, the corridor 
width is reduced. Figure 2.3-3 shows a plot of the control 
corridor width versus the L/D ratio. This figure represents the 
maximum corridor width available with a perfect navigation and 
guidance/control system, and with no orbital inclination change. 

The minimum control corridor width required, resulting from 
navigation (entry state) error and aerodynamic uncertainties, is 
summarized in Figure 2.3-4. The detailed analysis for each of the 
variables can be found in the sensitivity study (Sections 5.2 and 
5.5) . Combining the effect from the listed errors and 
uncertainties yields a minimum corridor width of 4.99 nautical 
miles. A 35% margin results in a control corridor width 
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requirement of' 6.75 nautical miles. 


Figure 2.3-5 shows the impact of the L/D ratio on the maximum 
plane change capability. It indicates that the higher the lift the 
more turning capability is available. However, an orbit plane 
change reduces the available control corridor width because the 
vehicle must spend time to correct direction as well as speed. The 
effect of the plane change can be seen in Figure 2.3-6. 

2.4 Guidance System (HYPAS) 

The HYPAS (Hybrid Predictive Aerobraking Scheme) guidance 
algorithm was obtained from MSFC and implemented in the 6D0F and 
3D0F simulations. The scheme consists of two coupled algorithms 
for aeroassisted orbit maneuvers. One algorithm operates to 
control the out-of-plane motion in order to achieve the desired 
orbit plane. The second algorithm operates in the vertical plane 
by targeting the exit conditions to achieve the desired orbit 
apogee. The vertical plane algorithm has two phases. The first 
phase is an equlibrium glide phase in which the guidance 
parameters can be adjusted to control the minimum altitude, 
maximum g-load, and heating rates. The exit phase (second phase) 
guides to the desired apogee altitude and attempts to minimize the 
AV required for insertion into the final low earth orbit (LEO) . 

The algorithm output is a bank angle command. The angle-of-attack 
for the lifting brake is constant and the bank angle is used to 
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control the direction of the lift vector. A negative bank angle 
command will pull the vehicle to the left and a positive bank 
angle command will pull the vehicle to the right (velocity axes) . 
If the magnitude of the command is less than 90 degrees, the lift 
vector will pull the vehicle away from the Earth. If the 
magnitude of the command is greater than 90 degrees the vehicle 
will be pulled toward the Earth. The four quadrants are shown in 
Figure 2.4-1. 

2.5 Control Technique 

A reaction control system was developed to manuever the vehicle in 
response to the guidance commands. A block diagram of the AOTV 
classical control system and attitude control blocks is shown in 
Figure 2.5-1. The reaction control systems reside in the pitch, 
yaw, and roll blocks. The inputs to the RCS are the attitude 
errors and the outputs are torques about the corresponding axes. 
The following two sections describe the choice of the control 
effector and the operation of the control system. 

2.5.1 Control Effectors 

Two types of control effectors were considered for the AOTV, 1) 
reaction control jets (thrusters), and 2) aerodynamic trim 
surfaces. Thrusters are required for control in the space 
environment, while the aerodynamic trim tab could be used in the 
aeropass. The added complexity of the trim tab (control surface, 
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actuator, power source) and the lack of a requirement for precise 
angle-of-attack control reduced the need for the trim surface. 
The guidance algorithm does not require control of the 
angle-of-attack because it assumes that the vehicle is stable and 
operating at the aerodynamically trimmed angle-of-attack. It is 
not feasible to trim the vehicle a different angle-of-attack with 
the reaction control jets in the atmosphere due to excessive fuel 
use. 

The thrusters were located such that they do not impinge upon the 
vehicle' structure and are shielded from the heating wake. The 
locations are at 15.7 feet fore and aft of the center of gravity, 
and 5.0 feet on each side of the centerline. Figure 2. 5. 1-1 shows 
the vehicle and the thruster locations. Only the yaw thrusters 
create a pure moment on the vehicle; the pitch and roll thrusters 
also produce a force along the positive body Z axis. At the 
trimmed angle-of-attack, this produces a small AV effect. 

2.5.2 Control System 

Two RCS control techniques were developed for the AOTV. Both 
methods were designed to accommodate bank and angle-of-attack 
commands from the guidance algorithm. The first control design 
uses classical linear control theory with proportional thrusters. 
Recent advances in thruster technology have made the use of 
throttleable thrusters feasible. The control system uses 
attitude and derived rate feedback in a proportional/derivative 
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control architecture. The gains for each axis were chosen to give 
a 0.2 Hz response and a damping ratio of 0.8. 

The second design uses conventional on-off thrusters in a 
phase-plane control algorithm. • The phase-plane uses parabolic 
switching curves derived from the maximum acceleration capability 
and a rectangular dead zone to mimimize minimum-bit limit cycling. 
The parabolic switching curves give minimum time respone to an 
error signal and no chattering along the switching line. Figure 
2. 5. 2-1 is a block diagram of the two control system layouts. The 
pitch, yaw, and roll control law blocks in Figure 2.5-1 can be 
either of the above mentioned control techniques. Figure 2. 5. 2-2 
lists the parameters for each axis for both control techniques. 
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3 . 0 DESIGN CRITERIA 


A set of design criteria, such as the stabilty margins and 
response times, are usually specified for the control design. The 
6D0F simulation with atmosphere and vehicle/hardware uncertainties 
is then used to evaluate the performance and robustness of the. 
design. This section details the design criteria that were used to 
develop the AOTV control system. 

3.1 Fuel Use And Control Response Requirements 

The criteria for the design of the AOTV autopilot were to create a 
control design that would minimize the fuel burned during the 
atmospheric maneuvers and also to closely follow the guidance 
commands to minimize the fuel needed for insertion into the 
desired orbit. These tend to be opposing criteria since the first 
objective requires minimal thruster use and the second objective 
requires significant thruster use to follow rapidly changing 
large angle guidance commands. 

The autopilot must have sufficient bandwidth and accuracy to 
insure that the total system achieves a desirable post-aeropass 
velocity and trajectory. An indication of the robustness can be 
identified in the gain and phase margin of the system. These 
margin are calculated by opening each loop of the autopilot and 
analyzing the stability margin. The design goal for each autopilot 
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loop, opened at the actuator command, is +6 db gain margin and +30 
degrees phase margin. The overall system robustness is verified 
with the sensitivity analysis, wich is measured in terms of fuel 
consumption and DV. 

3 . 2 Atmosphere Characteristics 

Variations in the atmosphere directly affect the AOTV mission. The 
control system must be able to handle variations in density by as 
much as twenty percent. The system must also accommodate winds of 
up to 400 feet per second. Figure 3. 2.-1 shows sample air density 
and wind profiles generated for each season using the GRAM 
atmosphere model (Global Reference Atmosphere Model) . 

3.3 Design Uncertainties 

The control system must be able to handle uncertainties in the 
mass distribution, aerodynamic models, and sensor outputs. The 
following section lists an estimate of the vehicle and hardware 
uncertainties. The uncertainty values are the result of past 
experiences on other projects and associated hardware. These 
values are used in the sensitivity analysis in Section 5. 

3.3.1 Mass Properties 

In the design phase, the uncertainty in mass properties is large. 
As the project evolves to the final design and the payload and 
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mass properties are more accurately determined, the uncertainty 
value approaches 2-5%. Each payload will alter the c.g. location 
of the vehicle, which will in turn alter the trim position. The 
control system must remain stable in the c.g. -shifted condition. 

3.3.2 Aerodynamics 

Table 3.1 lists estimates of the aerodynamic uncertainties. These 
tolerances apply to a vehicle which has been aerodynamically wind 
tunnel tested. The values are associtated with the wind tunnel 
test accuracy. The uncertainties are small (±5%) for force data, 
with somewhat larger uncertainty for the moment data. The 
tolerances are large for the moment cross-coupling term since the 
nominal values are usually small. 

3.3.3 Strapdown Platform Uncertainties 

Table 3.2 lists the tolerances associated with a strapdown 
platform. Values shown are specification values received from 
strapdown platform vendors. These values are for state-of-the-art 
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4.0 DESIGN SYNTHESIS 


This section presents the results of the trade studies used to 
determine the final control system design. The results were used 
to modify and fine tune the control system and the guidance 
algorithm to achieve better mission performance. 

4.1 Thruster Sizing 

Thruster sizes of 100 lbs, 250 lbs, 350 lbs, and 500 lbs were 
traded to determine the required size. Initially, all axes (roll, 
pitch, and yaw) used the same size of thruster. In the final 
design, smaller thrusters were used in the pitch axis because the 
bank commands from the guidance algorithm involve only the roll 
and yaw thrusters (for trimmed angle-of-attack) . 

Rate limits of 15 deg/sec, 17 deg/sec, 20 deg/sec, and 30 deg/sec 
were used in the autopilots. The rate limits were used to smooth 
the bank command and eliminate excessive fuel use in the 
thrusters. Figure 4.1-1 is a block diagram of the digital 
autopilot design. 

The effects of the RCS thruster size and the autopilot rate limit 
were examined using the total fuel and total AV as measures of the 
sensitivity. Figure 4.1-2 shows the fuel use versus thruster size 
for all four rate limit configurations. The 15 deg/sec case 
indicates the lowest fuel use but in fact the lag was so extreme 
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with respect to the guidance algorithm that the autopilot 
effectively did nothing. This is reflected in Figure 4.1-3 for the 
15 deg/sec curve which shows large AV requirements for any 
thruster size. Larger thrusters combined with the larger rate 
limits performed better at the expense of' more fuel burned. A 
somewhat "optimal" design can be chosen by trading fuel used for 
required AV. The 20 deg/sec rate limit gives the lowest AV for 
the smaller thruster sizes and the median response for the fuel 
use. A 250 lb thruster gives a Av of around 500 ft/s and uses 
around 250 lbs of fuel. This design point is shown on Figures 

4.1- 2 and 4.1-3 by X's. 

The operation of the guidance algorithm is the reason for the 
seemingly large thrusters and large amount of fuel use. The bank 
commands from the algorithm require large (-165 deg, +165 deg), 
fast rotational maneuvers of a vehicle with large rotary inertia 
(74500 slug-ft 2 ) . The smaller thrusters are just capable of 
following the guidance close enough to reduce the exit condition 
errors . 

Figues 4.1-4, 4.1-5, 4.1-6, and 4.1-7 are time history responses 
of the 6DOF simulation using the design point with 250 lb 
thrusters and a 20 deg/sec rate limit in the autopilots. Figures 

4.1- 4 and 4.1-5 show the altitude versus velocity and altitude 
versus time trajectories. Figure 4.1-6 is the bank command and 
actual bank angle. Notice that the actual bank angle lags behind 
the commanded bank angle. Figure 4.1-7 is the response of one of 
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the yaw thrusters . 


The final design uses 250 lb thrusters for the yaw axis and 100 lb 
thrusters for the pitch axis. The rate limits were established at 
20 deg/sec. 

4 . 2 Control System Trades 

Two types of thruster/autopilot systems were included in the six 
degree-of-freedom simulation, a proportional thruster/autopilot 
system and a phase-plane/autopilot system. A proportional system 
is more effective for studying the vehicle sensitivities and the 
guidance algorithm because of its response characteristics. The 
proportional system contributes to the dynamic response in a 
predicable way which makes the other system contributions more 
distinguishable . 

There is an inherent problem created when using a 
phase-plane/autopilot system and a rapidly changing command source 
such as HYPAS . Phase-plane systems are usually designed to 
accommodate a step position command, not follow a series of 
commands which have faster time constants than the vehicle 
response time. One way to accommodate these types of commands in- 
a phase-plane system is to use large thrusters to lower the 
vehicle response time. Without the large thrusters in the roll 
and yaw axes, the vehicle was not able to follow the guidance well 
enough to sufficiently reduce the velocity errors at exit. 
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The phase-plane was designed using parabolic switching curves and 
a rectangular dead zone for position and rate errors. This 
phase-plane design always used more fuel than the proportional 
autopilot for the same flight scenario. Figures 4.2-1, 4.2-2, 
4.2-3, and 4.2-4 show the. phase-plane verus proportional autopilot 
fuel trade. Given a specific vehicle design, weight, thruster 
size, etc., a more elaborate phase-plane could be developed which 
would approach the proportional autopilot performance. 

4 . 3 Adaptive Angle of Attack Control 

An adaptive control scheme was designed to "hunt” for the 
aerodynamic trim condition in the pitch axis as the 
center-of-gravity was moved with different payloads. The 
integrator adjusts the commanded angle-of-attack until the vehicle 
is trimmed. This would reduce the amount of fuel used by 
eliminating the need for the pitch thrusters to counter the moment 
caused by the non-trimmed flight condition. 

The technique worked well for some C.G. locations and not for 
others. In fact, for some C.G. locations the vehicle did not even 
exit the atmosphere (mission failure) . Figure 4.3-1 can be used 
to explain the effect of the trim integrator on mission 
performance. The figure shows a plot of the coefficient of lift 
as a function of angle-of-attack and plots of the trim values of 
the vehicle for x and z variations in the C.G. location. Notice 
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that as the x or z location of the C.G. was increased (positive 
direction) the value of Ci corresponding to the new trim condition 
was decreased. Thus at the new trim condition there is 
significantly decreased lift available to modulate the in- and 
out-of-plane errors in the trajectory. As the C.G. locations were 
moved in the negative direction the corresponding values of the 
lift coefficient increased. 

The gains in the guidance algorithm were also modified to be 
semi-adaptively adjusted as a function of the estimated dynamic 
pressure. The equations were obtained from documentation on the 
derivation of the HYPAS guidance scheme [1] . These equations 
maintain an overdamped (zeta = 1.5) 0.05 hz response in the 
guidance logic. Other values for the damping and response 
frequency were tried, but all decreased the performance of the 
system. The damping ratio was decreased to zeta = 1.0 and the 
response frequency increased to 0.1 hz with mixed combinations of 
the two parameters in between. The addition of the equations to 
the simulation made only a small improvement in the system 
response. This is to be expected because most of the effective 
control authority occurs over a very small portion of the 
trajectory where the dynamic pressure is relatively constant. 

4.4 Guidance Law Verification and Sensitivity 

Initial comparisons of the 6D0F and 3DOF simulations did not match 
well. Specifically, the required AV' s were much larger for the 
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6D0F simulation. After much tuning and verification of the 
algorithm, the response was improved but under no circumstances 
was able to match the performance of the 3D0F version. The basic 
difference between the two is that the rotary inertia in the 6D0F 
simulation prevents an exact following of the bank commands. 
Larger thruster sizes obviously improve the response, but at the 
expense of more fuel use. 

A rate limiting first-order filter set at 0.2 Hz with a 30 deg/sec 
rate limit was used on the bank command output to smooth the 
commands. Without the filter the commands drove the RCS 
autopilots very hard and used excessive fuel. This situation 
occurs when the guidance alogrithm is not on a hard limit, and 
issues a ratchety command sequence. 

The best performance was obtained by increasing the gain on the 
vertical plane correction during the entry phase, then decreasing 
the gain in the exit phase. This effectively creates the 
situation where the vertical plane components are corrected in the 
entry phase and the out-of-plane components (inclination) are 
corrected in the exit phase. This method decreased the required 
AV after exit by several hundred feet per second depending on the 
autopilot' and guidance parameters selected. 

The following are the parameters used in the guidance loop for the 
preliminary sensitivity studies: 
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GHDOT =*.40 


GHDOTX =*.25 
VSAT =26700 ft 
VTRIG =29000 ft 
AMXER =.63 deg 
GN17b =7.6 


These parameters were used as the design point for the guidance 
algorithm for all sensitivity studies. 

The HYPAS guidance algorithm was further modified to save fuel by 
setting the command to zero during the exit phase when the dynamic 
pressure dropped below .5 lbs/ft 2 . Aero maneuvers beyond this 
point did little to reduce the exit errors. 


40 



Rate Limited Filter 



4 # 


<T) CO 
C ® 
® <? 
CO <0 


Figure 4.1~1 Digitdl Autopilot Filter 












43 





HEWLETT-PACKARD 92701027 



HEWLETT-PACKARD 9270-1027 



45 


ORIGINAL PAGE IS 

OF POOR QUALITY 










(63p) ^UHQ 


46 

ORIGINAL PAGE IS 

OF. POOR QUALITY 



HEWl.liTT-PACKARO 9270 1027 










number 

REV LTR 



01-40138001 REV. 11/84 


ORIGINAL PAGE IS 
OF POOR QUALITY 




ORIGINAL PAGE IS 
OF POOR QUALITY 


NUMBER 
REV LTR 



01-4013-8001 R6V. 11/84 


SHEET 

49 



&£7&JSV£Z 


NUMBER 
REV LTR 



01-4013-8001 REV. 11/84 


SHEET 

50 


ORIGINAL PAGE; !S 
OF POOR QUALITY 


*i.70 



01 4013-8001 B6V. 11/84 


SHEET 

51 


Figure 


O * 
n % 




CO ° 

i 

oo 

1 

l 


*5J(7 ~ **'^70 





□ i4oi 38001 nev 11 m* 


SHEET 

52 


OC ^ Oes. 








5.0 SENSITIVITY STUDIES 


The 6DOF simulation, including the guidance and control 
algorithms, was used to evaluate the sensitivity of the AOTV 
mission to uncertainties in six major categories: 1) mission, 2) 
entry state, 3) mass properties, 4) atmosphere and winds, 5) 
aerodynamics, and 6) sensors. The sensitivity analysis for each 
of the six areas is explained below. Figure 5.0-1 is a summary of 
the sensitivity analysis. 

Sensitivity was measured in three ways. The first measure was the 
amount of fuel burned during the aeropass. The second measure was 
the size of the AV burns necessary to establish the desired orbit. 
The third sensitivity measure is the combination of the previous 
two, yielding a total sensitivity measure. 

The three AV burns are shown in Figure 5.0-2. The first burn 
raises perigee, the second burn circularizes the orbit, and the 
final burn corrects the inclination. Figure 5.0-3 shows the 
subroutines and logic used to calculate the required burns. 

The fuel required for the AV burns can be calulated given the mass 
properties of the OTV the specific impulse of the thruster fuel, 
and the required velocity change. A AV burn corresponds to 
pounds of fuel through the relationship lb f = (m/I sp ) *AV, where lb f 
is pounds of fuel, m is mass of the vehicle, I sp is the specific 
impulse of the fuel, and AV is the required velocity change. 
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Using a specific impulse of 250 seconds, and a mass of 22000 slugs 
we get the formula lb f =2.73*AV. Thus each foot per second of 
velocity change requires 2.73 pounds of fuel. An effective 
guidance/control system will try to minimize the AV burns and 
thereby save fuel. 

The first sensitivity measure was the amount of RCS fuel used 
during maneuvers in the atmosphere. The amount of fuel burned in 
the aeropass can be converted to a AV as explained in the previous 
paragraph. Then the AV can be compared to the velocity change of 
the vehicle due to the atmospheric braking effects. The ratio of 
the AV due to thrusters to the AV due to atmospheric effects can 
be used to measure the effectiveness of the aerobraking scheme. 
The smaller the ratio the better. For example, a typical mission 
burns around 225 pounds of fuel with a velocity drop of 9500 feet 
per second. The fuel translates into a AV of 82.4 feet per second 
and the ratio becomes 0.0086, meaning that the aerobraking 
maneuver was an effective method of removing velocity. 

Parameter variations were accomplished by either adding 
incremental changes to the variable of interest, or multiplying 
the variable by a scale factor. 

Two FORTRAN subroutines are used to calculate the trajectory 
orbital parameters and the AV burns required to establish the 
desired orbit. These subroutines and an input/output map are 
listed in the appendix. Subroutine ORBIT may be called at any 
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time during the simulation to determine the orbital condition. It 
must be called when the vehicle reaches the exit condition and 
previous to the call to subroutine DVEL. Subroutine DVEL will 
calculate the three trajectory correcting velocity burns which 
will be required to position the vehicle in the desired orbit. 
These calculations are based on Hohmann Transfer Orbits followed 
by a plane correction. The subroutine outputs are used as one 
sensitivity measure, along with the total fuel used during the 
atmospheric phase. 

5 . 1 Mission Sensitivities 

Two mission profiles were considered. The nominal mission enters 
the aeropass with no inclination error and attemps to maintain the 
same inclination in the presence of the rotating earth. The 
second trajectory enters the aeropass with a 1.5 degree 
inclination offset and attempts to change the orbital plane to the 
desired value of 28.5 degrees. The closed and open-loop missions 
are initiated with the same state vector (total velocity, 
altitude, heading, flight path angle, inclination) . The nominal 
mission begins with no inclination error. The nominal trajectory 
results in large rotational motions to correct small inclination 
errors. This action burns the majority of the fuel used during 
the aerobraking phase. 

The second trajectory was generated using the 3D0F simulation, 
with a 1.5 degree initial inclination error. The intent was to 
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find a trajectory which does not start with zero inclination error 
and examine the system performance. The guidance tries to 
eliminate the inclination error, which requires the vehicle to 
remain on one side of the bank plane for the majority of the 
mission. Fuel use was reduced from 225 pounds for the nominal 
system to 161 pounds for the out-of-plane trajectory. Figure 

5.1- 1 is a plot of the bank angle and bank angle command for the 
out-of-plane trajectory. Notice that the bank angle and command 
are one-sided. Figure 5.1-2 shows the orbital inclination during 
the aerobraking phase. The desired inclination (28.5 degrees) is 
reached during the late phase of the mission. 

All remaining sensitivity studies were conducted with the nominal 
trajectory (no initial out-of-plane error) . 

5.2 Entry State Sensitivities 

The mission sensitivity to total velocity, altitude, inclination, 
and flight path angle errors are shown in Figures 5.2-1, 5.2-2, 

5.2- 3, 5.2-4, respectively. The left-hand axis and the dashed 
line on the plot correspond to AV measurements. The right-hand 
axis and the solid lines correspond to the RCS fuel used during 
the aerobraking phase. Each of the parameters exhibit similar 
sensitivity patterns, a relatively flat zone in the middle with 
rapid divergence when the parameter is too far off nominal. The 
sensitivity bounds are approximately ±6% off nominal for the 
velocity, ±8% for altitude, ±.4% for inclination, and ±9% in 
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flight path angle. By far the most sensitive parameter is the 
orbital inclination. 

5.3 Mass Property Sensitivities 

The mission sensitivities to uncertainty in vehicle weight, 
inertia, and x and z c.g. location are shown in Figures 5. 3.-1, 
5.3-2, 5.3-3, and 5.3-4. 

Variations of the weight by as much as ±27% resulted in changes in 
fuel use by at most ±23% and in the Av measurement by at most 
±17%. in summary, while the mission performance is affected by 
the vehicle weight, it is not nearly as sensitive to weight as it 
is to the entry state parameters. 

Increasing the inertia resulted in a proportional increase in the 
fuel. The AV curve should be viewed in light of the scales. A 
±40% variation in the inertia resulted in. a ±2.4% variation in the 
AV measurement. In summary, the mission is not highly sensitive 
to inertia and the effect is predictable. 

Figures 5.3-3 and 5.3-4 are plots of the sensitivities to x and z 
variations in the c.g. location and will be discussed together. 
The two solid lines in each plot are the fuel used with and 
without the integrator logic. The two dashed lines are the AV 
measurements with and without the integrator measurements. From 
the figure we can see that the mission is least sensitive to 
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variations in the x direction without the integrator. Addition of 
the integrator increases the performance if the c.g. is moved in 
the negative direction and causes the system to fail if the c.g. 
is moved too far in the positive direction. Too much variation of 
the c.g. in the z direction (either direction) causes the system 
to fail without the integrator. The addition of the integrator 
increases the performance if the c.g. is moved in the negative 
direction and quickly causes the system to fail if the c.g. is 
moved in the positive z direction. This is easily explained by 
looking at the aerodynamic data and Figure 4.3-1. If the c.g. is 
moved in the negative direction, the trimmed angle of attack 
decreases, which results in an increased lift capability. 

5.4 Atmospheric Sensitivities 

The atmospheric sensitivity studies included the effect of density 
variation and winds. Sixteen random atmosphere models were 
generated using the GRAM (Global Reference Atmospheric Model) 
atmosphere program. The models include random variations in 
density and wind profile. Figure 5.4-1 is a plot of the density, 
north-south, and east-west winds as a function of altitude. The 
sixteen atmospheres were chosen as four worst cases from each 
season - summer, fall, winter, and spring. The four cases for 
each season are maximum dynamic pressure (MQ) , worst winds (WW) , 
most northerly (MN) , and most southerly (MS) . The nomenclature 
used to identify the winds is a four letter name with the type of 
model first and the season as the second two letters. For example 
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a most southerly fall profile would be MSFA. Summer is (SU) , 
winter is (WI) , and spring is (SP) . Figure 5.4-2 shows the fuel 
and AV numbers for the nominal trajectory using the sixteen 
atmospheres, with and without winds. No general trends are 
discernible because the atmosphere models are random in nature. 
However, the vehicle remains stable and completes the mission 
profile for all of the different atmosphere models. 

5 . 5 Aero Coefficient Sensitivities 

The mission sensitivity to uncertainty in the aerodynamic 
coefficients was determined by multiplying the aero coefficients 
by a scale factor. The nominal scale factor value was one. 
Sensitivity studies were conducted on the force (C x , C y , C z ) , 
static stability (Cm a , Cn B ) , and yaw/roll coupling (C1 B ) aero 
coefficients. 

Figure 5.5-1 shows the sensitivities for the three force coeffi- 
cients. The horizontal axis has units of percent uncertainty from 
the nominal condition. The combined sensitivity measure clearly 
shows that the effect of off-nominal C y is negligible, as is 
negative variation of the C z coefficient. Positive variations in 
the C z coefficient are more pronounced but do not cause the 
control /guidance systems to fail. The C x coefficient has a stable 
sensitivity range of -15% to +5% off the nominal condition. 

Figure 5.5-2 is a plot of the Cn B , C1 B , and Cm a sensitivities. 
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Notice that the C1 B and Cn B curves are overlayed. The vehicle is 
insensitive to these parameters . The two remaining sets of curves 
show the Cm a sensitivities with and without the pitch axis 
integrator. The integrator has a stabilizing effect for positive 
variations of the coefficient. For negative variations the 
integrator quickly causes the system to diverge. This is due to 
the decreased control authority at the reduced angle of attack. 
Figure 5.5-3 is a plot of the combined sensitivity measure for 
these parameters. The effect of the integrator can be seen very 
clearly. 

5.6 Sensor Sensitivities 

The mission sensitivity to noise and bias in the strap-down sensor 
package was evaluated by including a detailed model of a rate gyro 
package in the 6D0F simulation. 

Figure 5.6-1 shows the model for rate gyro errors. Five error 
sources have been modeled: misalignment of the rate gyro package, 
bias and scale factor error, additive noise in each channel, and 
the gyro cross-axis effects. The axes of the gyros in the package 
were assumed perfectly orthogonal, with the alignment error 
defined from the package roll, pitch, and yaw axes to the true 
body axis. Table 5-1 lists the expected 3a statistic value for 
each parameter. The listed values apply to missile systems. For 
space hardware, the requirements are more stringent because 
mission durations are usually much longer. 
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Figure 5.6-2 shows the sensitivity for the gyro. parameters except 
for noise. The plots show that the system is relatively 
insensitive to errors in the rate gyro. The amount of fuel burned 
due to the errors is small (less than 5 percent) . The rate gyro 
errors are principally manifested in the increase of AV 
measurement. The increase is comparatively small (less than 13 
percent) . 

The system is also insensitive to noise in the rate gyro channel, 
as the dynamic response of the vehicle acts as a low pass filter 
in removing most of its effect. The effect of gyro noise can be 
seen, however, in the increased thruster activity level, which 
translates to an increase in fuel burned during maneuver phase. 
The fuel burned due to 3a noise level in the roll, pitch and yaw 
rate gyro .channels are 269, 255 and 254 lbs respectively, compared 
to 244 lbs of fuel for a nominal flight. Again, the increase in 
fuel burned is small. 

In conclusion, the mission is most sensitive to parameters that 
directly affect the guidance alogorithm, such as inclination 
error. The control systems developed for the AOTV were stable for 
all parameter variations. The performance decreases were 
principally due to the minimal control authority which was only 
significant during the one to two minutes of highest dynamic 
pressure . 
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The HYPAS guidance algorithm is listed in the appendix 
reference and comparison. 
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Subroutine DELY 


inputs : 
outputs : 

Subroutine 
inputs : 

outputs : 


RP1 — > 
RA1 — > 
RCIRC — > 
RAAN — > 
RINC — > 


radius at perigee (feet) 

radius at apogee (feet) 

radius circular orbit (feet) 

right accension of acending node (deg) 

orbit inclination (deg) (28.5) 


DV1 — > 
DV2 — > 
DV3 — > 
DVEL — > 


.burn to correct perigee 
burn to correct inclination 
burn to circularize 
total delta velocity 


QRBI1 


Reci x Y z — > ECI position vector (feet) 
Veci,/,/- — > ECI velocity vector (feet /sec) 


RPER — > 
RAPO — > 
.RCIRC — > 


radius at perigee (feet) 
radius at apogee (feet) 
radius of circular orbit (feet) 


Figure 5.0-3 cont. 
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SUBROUTINE ORBIT 
C 

SAOO V : AOTVCM 
C 

C MOMENTUM VECTOR AND MAGNITUDE 
HX»RECIY»VECIZ-VECIY«RECIZ 
HY»RECIZ*VECIX-VECIZ*RECIX 
HZ*RECIX*VECIY-VECIX*RECIY 
HMAG«SORT( HX*»2*KY**2+HZ»*2 ) 

C NODAL VECTOR AND MAGNITUDE 
RNX--HY 
RNY-HX 
RNZ»0.0 

RNMAG-SORTI RNX»»2*RNY**2*RNZ**2 ) 

C ECCENTRICITY VECTOR AND MAGNITUDE 

V$QRD»VEC I X**2*VEC I Y**2+VEC IZ**2 
RDOTV»REC I X*VEC I X+REC I Y*VEC I Y+REC IZ*VECIZ 
TEMPZ»VSQRD-RMU/( -Z? ) 

EX»1 . / RMU* ( T£M?Z*REC I X-RDOTV*VEC I X ) 

EY-1 . /RMU"< TEMPZ*REC I Y-RDOTV*VEC I Y ) 

E2»I ./RMU»( TEMPZ*REC:Z-P.DOTV«VECIZ> 

EMAG«SQRT( EX«»2*E Y«»2*EZ*»2 > 

C APOAPSIS AND PERIAPSIS CALCULATIONS 
SLR»HMAG*»2/RMU 

I P ( HMAG .NE. 0.0) R I NC»ACOS< MZ/HMAG ) “Z57P 3 
I F ( RNMAG .NE. 0.0) RLAN-ACOSt RNX/RNMAG )»ZS7P3 
RPER-SLR/( I .0*EMAG 
RA?C*SLR/< 1 . 0-EMAG > 

RPEP.N«< RPE R-R EARTH .'/ S076 . I 155 
RAPOIW RAP0-REAR7H J/5C7S .1155 
RCIRC*REARTH*APOTG*6076. 1 155 
RETURN 
END 


Figure 5.0-3 cont . 
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SUBROUTINE OELVt RP 1 . RA 1 . RC I RC . RI NC . RAAN . DV1 , DVZ. 0V3 . OVEL ) 

XMU ! GRAVITATIONAL PARAMETER IN FT»»3/SEC 

DATA XMU / 1 . 4076E+ 1 6/ 

DATA 257P3 /57. 29578/ 

EXIT ORBIT VELOCITY AT APOGEE: 

VAl-SORT(2.-XMU"< I/RA1-1/<RP1+RA1 > ) ) 

TRANSFER 0R8IT VELOCITY AT APOGEE: 

VA2*SQRTC2.*XMU*( 1 / RA1 - 1 /< RA1 + RC I RC ) ) ) 

DELTA V REQUIRED FOR EXIT TO TRANSFER ORBIT CHANGE: 

DV1-ABS( VA2-VAI ) 

CIRCULAR ORBIT VELOCITY.: 

VCIRC*SQR7( XMU/RCIRC ) 

TRANSFER ORBIT VELOCITY AT PERIGEE 

VP2-SQRT< 2.*XMU*( l/RCIRC-l/< RAl+RCIRC) ) ) 

DELTA V REQU I RED FOR TRANSFER TO CIRCULAR ORBIT CHANGE: 

DV2»A8S( VCIRC-VP2 ) 

CALC PLANE CHANGE DELTA V {TARGET ORBIT INCLINATION AND RIGHT 
ACCENS I ON OF ACENDING NODE ARE HARD CODED AT 28.5 DEG AND 0 
RESPECTIVELY 

RINCR-RINC/ZS7P3 
RINCRF-28.5/Z57P3 
RAANR-RAAN/Z57P3 
IF {ABS(RAAN) .LT. I.E-8) THEN 
PL NCHG»AB S (RINCR-RI NCRF > 

ELSE 

PCCOS«COS<RI NCR) "COS <RINCRF)*SIN(RINCR)*SIN(RI NCRF )*COS( -RAANR ) 
PLNCHG-ACOS(PCCOS) 

END IF 

DV3*4BS( 2 . *VC IRC*S !N<PLNCHG/2. ) ) 

TOTAL DELTA V 

DVEL*DV1+DV2»0V3 

RETURN 

END 


Figure 5.0-3 cont 
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6.0 ADVANCED CONTROL SYSTEM CONCEPTS 

Advanced control concepts are used to obtain improved system 
performance and robustness, relative to established control 
techniques and designs. These advanced concepts may include an 
adaptive element which is used to identify unknown or uncertain 
parameters in the system. For example, an adaptive algorithm 
might fine tune a control system by identifying the actual drag 
coefficient, as opposed to using the value generated by the 
aerodynamic model . 

Another group of advanced control techniques includes the 
so-called "modern methods, " which address the multiple-input, 
multiple output control problem. The thrust of these techniques 
is generally to minimize some form of a cost function, which 
varies for each technique. A well known example is the quadradic 
cost function of the Linear Quadradic (LQ) technique, which may be 
configured to minimize time, fuel, distance, or other important 
parameters . 

♦ 

The object of the aerobraking scheme is to maximize the energy 
"bleed-off" during the aeropass while minimizing the exit window 
error such that minimum fuel is required for orbit correction 
burns. An additional objective is to minimize the fuel burned 
during the aeropass. The first objective translates into a 
guidance requirement, i.e., use of an algorithm that will 
"bleed-off" energy and minimize exit window errors. The second 
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objective can be addressed using classical or modern control 
methods . 

The following paragraphs describe the design and implementation of 
a semi-adaptive pitch axis controller. A new adaptive scheme 
which control combines aspects of each the techniques to achieve 
improved performance, robustness, and efficiency is presented in 
Section 6.2. 

6.1 Semi-adaptive Pitch Axis Control System 

The pitch axis rotational dynamics can be characterized with the 
second order transfer function 

£= Mg 

8 <s 2 - M a ) 

where 9 is the angular position and Sis the commanded angular 
acceleration. M a and M s are the moment sensitivities, which 
for this vehicle reduce to 

M a = Cm a *Q*S*D*57.3 = -K 

lyy Iyy 

m 6 = 57 .3 

Iyy 
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where 57.3 is used to keep the units consistent (degrees) . K may 
be viewed as the system "stiffness" which varies with dynamic 
pressure. 

The natural frequency is given by the square root of (defined 
above) . Over the course of the AOTV mission the natural frequency 
ranges from 0.0 to 2.1 radians per second, a function of the 
dynamic pressure. An autopilot designed to keep the pitch axis 
well damped at low Q will have a decreased damping ratio as Q 
increases, see Figure 6.1-1. The semi-adaptive control scheme 
estimates dynamic pressure using the X or Z axis accelerometer 
measurement and the corresponding aero coefficient, and then 
adjusts the control gains and bandwidth to keep the vehicle 
damping ratio at .8 throughout the entire mission. At the large 
angle of attack (72 or 74 degrees) the X measurement is 
principally a lift measurement and the Z measurement is 
principally a drag measurement. Figure 6.1-2 is a block diagram 
of the dynamic pressure estimation logic and the semi-adaptive 
gain adjustment. The X axis acceleration measurement is depicted. 

The bandwidth of the semi-adaptive scheme is shown in Figure 6.1-3 
for the nominal mission. As the bandwidth increases the pitch 
axis thrusters are more active and more fuel is used. The 
semi-adaptive scheme uses 56 pounds more fuel than the nominal 
fixed gain controller (increase of 25 percent) . The semi-adaptive 
scheme also causes the required AV burns to increase from 521 feet 
per second to 530 feet per second (increase of 1.7 percent) . This 
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increase comes from the additional Z axis loadings by the 
thrusters. 

Figure 6.1-4 shows the pitch axis thruster firings for the 
adaptive and fixed gain autopilots.’ 

6.2 Adaptive Control Concepts 

The goal when using adaptive control techniques is to obtain 
acceptable levels of performance and stability with a changing or 
unknown plant. In the case of the AOTV vehicle, the parameters of 
most concern are the atmospheric variations and the c.g. shifts 
due to payload uncertainty. The atmospheric variations affect the 
guidance (HYPAS) and control, which both use an open-loop 
estimation scheme (see Figure 6.1-2) for dynamic pressure based on 
assumed knowledge of the aerodynamic coefficients. The c.g. 
shifts cause the trim position of the vehicle to change, requiring 
excessive thruster firings to maintain a specified trim angle. 
These two parameters represent the largest uncertainties and 
sensitivities in the system, and are therefore the most viable 
parameters to be identified in an adaptive control scheme. 

Model Reference Adaptive Control (MRAC) is most applicable to 
rmin i s t i c tracking problems, i.e. , when the desired output 
responses are known functions. For example, an advanced fighter 
can be made to handle like a training aircraft using 
MRAC. in this case, the desired response (trainer) is known 

90 


c- 1 



through on-board simulation, and the difference between the actual 
response (fighter) and the desired response can be calculated and 
used to influence the control, see Figure 6.2-1. The MRAC 
technique tries to minimize the error in the controlled variables, 
which for the example might be the normal acceleration, angle of 
attack, or pitch rate. Since there is no defined response function 
for the shaped-brake vehicle, this method is not recommended for 
the AOTV. 


Self Tuning Control (STC) modifies the control as the unknown or 
uncertain parameters in the plant are identified. The technique 
attempts to achieve improved control through minimization of a 
performance index. Example cost functions may include 
minimization of any function of the system parameters; for example 
minimum fuel use or minimum deviation from a desired set point. 
The STC methods usually include an identification algorithm such 
as the "least squares technique" which must identify the unknown 
system parameters in real time. These algorithms can be greatly 
simplified by keeping the parameters to be estimated to a minimum. 
Variations of this technique are recommended for the AOTV vehicle 
because useful cost functions such as minimum fuel are important 
to mission performance. 

The following is a proposed adaptive control design for the AOTV 
vehicle which uses a variation of the STC technique combined with 
a second non-adaptive controller and an optimization algorithm. 
Under this approach, a non-adaptive control system is designed 
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using classical or modern techniques (Linear Quadratic, 
H-Infinity, Eigenstructure, Optimal Projection, etc.) using 
expected values of the vehicle and mission parameters. This 
system is simulated and analyzed for stability. This is 
considered the base design. 

For each specific mission there may be new parameters/constraints 
such as payload requirements, atmospheric conditions, or maximum 
thruster size. The base design is then altered to accommodate the 
new parameters using the nonlinear optimization program QDES. 
QDES (Q-Design) is a computer aided design tool developed at 
Stanford by Boyd [2] . The QDES program reformulates the control 
problem using coprime factorization and searches for an optimal 
solution in the defined convex search space. If an optimal 
solution exists, the algorithm is guaranteed to converge to the 
global optimum.. The constraints which may be included in the 
optimization process are frequency domain type constraints such as 
phase and gain margins, bandwidth, and settling time. Using this 
program, the base design can be augmented in an optimal way to 
handle the new requirements. 

Up to this point no adaptive techniques have been incorporated. 
To this end, a Self Tuning Controller can be included in the above 
design for in-flight fine tuning of the control as better 
estimates of the true plant become available. The advantage of 
using the ST to fine tune the optimized base design over a purely 
adaptive control scheme is that the base design is already a 
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stable control system and the adaptive system does not need to 
identify the system to maintain stability. Thus the adaptive 
portion of the algorithm only fine tunes an already stable design 
for improved performance. 

The basis of the self tuning algorithm is similar to that of the 
QDES program, both being derived from a parameterization of the 
class of all stabilizing controllers as found in Vidyasagar [3]. 
Moore [4,5] has extended the parameterization technique to the 
class of all stabilizing, time-varying controllers for 
time-varying plants. Using the techniques of Moore, the stable 
transfer function matrix designed with the QDES program can be 
augmented with a time varying stable transfer matrix, which must 
be calculated in real-time. 

This advanced adaptive scheme has several advantages over a- purely 
adaptive control scheme. First, the knowledge base of past 
experience can be built into the base design. This provides a 
fail-safe level of control which can be used as a default if the 
adaptive algorithm fails to perform. Second, the control system 
is easily modified for different mission requirements without 
requiring extensive analysis for stability and performance. 
Third, the self tuning adaptive element attempts to tweak the 
control system to obtain optimal performance during the mission, 
but is not required to identify a large number of parameters and 
also guarantee stability. 
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The diagram in Figure 6.2-3 shows the structure of the proposed 
control system. Notice that Q and Q' can be zeroed out and the 
remaining control is the base feedback design closed through the X 
and Y _1 blocks (coprime factors) . With a nonzero Q and a zeroed 
out Q' the system has been augmented to handle the modified 
performance requirements. Q' is implemented in real-time and is 
the adaptive element in the system. This evolved control design is 
utilized to accommodate unknown system and mission parameters and 
achieve an acceptable level of performance. 
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Figure 6.1-2 cont. 
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Figure 6.6-2 Semiadaptive Control Based on Dynamic Pressure 










7.0 Conclusion 


Control technology has been developed in this contract for the 
AOTV mission which minimizes the fuel used for atmospheric 
maneuvering. The control techniques developed can be used on a 
number of AOTV type vehicles. The control algorithms can be 
tailored to the guidance law and vehicle characteristics by 
adjusting only the control gains and break frequencies. 

The control technique includes: 

1) Adaptive trim design for c.g. control; 

2) Semi-adaptive gain computer to minimize mission data loads; 

3) Proportional or on-off thruster selection; 

4) Control designs which minimize RCS fuel usage. 

Initial control studies have also been conducted on advanced 
adaptive control concepts. The concepts considered include 
model reference adaptive, linear quadratic, Q-design, and 
self-tuning adaptive. Results indicate that the classical 
design is sufficient for known vehicle configurations and 
payloads. Increased performance can be achieved with advanced 
control techniques over a large range of missions and payloads. 

The classical control techniques, have been verified through 
detailed sensitivity studies using 6DOF simulation. The 
simulation studies show the following: 

1) Small sensitivity (>20% tolerance) for all sensor effects. 
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CiP, CnP, Cy, and Cx aerodynamics, and mass properties 
(weight, inertia, X c.g. location) , 

2) Moderate sensitivity (10% to 20% tolerance) for Z c.g. 
location, atmosphere effects (winds, density) , and Cx and Cma 
aerodynamics; 

3) High sensitivity (<10% tolerance) for entry state parameters 
(velocity, flight path angle, altitude, and inclination) . 

The critical vehicle parameter in controlling the reentry state 
is the vehicle lift /drag ratio. The AOTV L/D is required to 
have a value larger than .15 to handle a 3o variation in the 
reentry state. Both vehicles analyzed, the 
and the AFE have L/D larger than .15. (L/D 
= .2, L/D AFE = .3) 


Boeing Shaped Brake 
Boeing Shaped Brake 
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8.0 RECOMMENDATIONS (FOLLOW-ON-STUDIES) 


The importance of the control and guidance techniques to the 
performance and robustness of the AOTV mission was studied in 
the Control Technology AOTV contract. New, advanced technology 
proposed for follow-on studies would benefit the AOTV by 
lowering recurring costs through autonomy and increase 
flexibility to varying mission requirements through adaptability 
and automated tuning. The autonomous design would make the 
system insensitive to payload variation and fuel load. Lower 
costs would be achieved by eliminating software changes per 
mission and minimizing mission data loads. The high 
reliability, low maintenance is achieved by using fault-tolerant 
processors and sensors. 

Our recommendations are to develop the technology for an 
autonomous three level modern control approach: 

1) Integral linear quadratic control design ( angle-of-attack, 
bank angle controller) 

2) Mission data load automation using Q-design method (tuned 
controller through simple mission data load) 

3) Adaptive fine-tuning based on parameter estimates generated 
in flight . 

This technology could be applied directly to the AFE 
configuration and verified on 6DOF simulation. Other guidance 
laws developed for the AFE by NASA Marshall Flight Center and 
Draper Labs can also be included in the 6DOF simulation. 
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APPENDIX 


A. 1 6D0F Simulation 

Figure A. 1-1 shows the hardware facilities used for the 6D0F 
simulation. The nonlinear equations of motion are hosted on the 
AD10. The autopilot and other slower processes are hosted on the 
Harris 1000. The AOTV simulation runs in real-time. 

Figure A. 1-2 is a block diagram of the entire 6DOF simulation. 
The dashed line indicates the separation of processes between the 
two computers. Each of the blocks in the diagram correspond to a 
FORTRAN subroutine in the simulation. 

A. 2 3DOf/6DOF Validation 


The 3DOF and 6DOF simulations are similar in that they use the 
same guidance algorithm. While there are many obvious 
differences, a few very important ones are discussed here. First, 
3D0F simulations neglect the effect of rotary inertia, an 
important quantity in the 6DOF simulation due to the rapid 
rotational commands from the guidance system. The 3DOF simulation 
uses a rate limiting scheme to approximate the rotational 
dynamics. We have found that a rate limit of 15 deg/sec best 
approximates the actual rotary dynamics for this vehicle. Second, 
the 3 DOF simulation uses a lift/drag aero model with constant 
coefficients while the 6DOF simulation derives the lift and drag 
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from the force and moment coefficients which are functions of 
angle of attack and sideslip angle. Third, the thrusters in the 
6D0F simulation create additional translational loads on the 
system that are not present in the 3DOF simulation. Fourth, the 
integration step size in the 3D0F is one second while the step 
size in the 6D0F is .001 seconds. Lastly, the scaling in the AD10 
computer used in the 6DOF simulation may have an effect on the 
results . 

Time history comparisons for the 3DOF and 6DOF simulations are 
shown in Figures A. 2-1 through A. 2-7. The figures show a good 
match in the simulations, verifying correct dynamic models and 
implementation of the guidance algorithm. Figure A. 2-7 is a plot 
of the bank commands. The difference in the two commands 
histories beginning at 160 seconds is due to the different 
calculated values of the inclination at that time. 

Table A. 2-8 shows the initial and final conditions for the nominal 
trajectory. The target window parameters are also included. 

The 3DOF and 6DOF simulations both use* the HYPAS guidance 
algorithm, which has two major modes of operation. During the 
entry phase the algorithm generates commands which attempt to 
maintain an equilibrium glide condition (lift = weight) . During 
the exit phase the algorithm generates bank commands based on the 
desired radius of apogee and the predicted exit condition. The 
difference between the actual orbit plane and the desired orbit 
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plane is used to determine the sign of the bank command. 
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SUBROUTINE GUID4 


SADO V : AOTVCM 
C 

C IMPLICIT DOUBLE PRECISION (A-H. 0-2) 

DIMENSION UYD(3),UYDS(3> 

REAL M , KRHOD . NRR 

C DOUBLE PRECISION M . MU . K 1 . KRHOD 

C DOUBLE PRECISION NRR 
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(AERO) 


SUBROUTINE INPUTS: 

NAVIGATION STATE: 

VR - EARTH RELATIVE VELOCITY (FT/S) 

BANK - CURRENT VEHICLE BANK ANGLE (DEG) 

ALF3IA - CURRENT VEHICLE ANGLE OF ATTACK (DEG) 

DRAGA « ACCLERATION DUE TO DRAG (FT/S~2> 

FP ACC * ACCELERATION VECTOR IN FLIGHT PATH COORDINATES 
( FT/S~2 ) 

X Y2I ■ ECI POSITION VECTOR (FT) 

XYZID * ECI VELOCITY VECTOR (FT/S) 

DESIRED END CONDITIONS: 

HEXIT - DESIRED EXIT ALTITUDE (FT) 

APOTG - DESIRED EXIT APOGEE ALTITIUDE (NM) 

DXY2I - DESIRED EXIT ECI POSITION VECTOR (FT) 

DXY2ID » DESIRED EXIT ECI VELOCITY VECTOR (FT) 

PHYSICAL CONSTANTS: 

MU - GRAVITATIONAL CONSTANTS <FT^3/S~2) 

RE ■ EARTH RAOIUS (FT) 

FPNM * FEET PER NAUTICAL MILE CONVERSION FACTOR 
RADCON « RADIANS PER DEGREE CONVERSION FACTOR 
VEHICLE PARAMETERS: 

CL ■ COEFFICIENT OF LIFT 

CD - COEFFICIENT OF DRAG 

VTZ - VEHICLE WEIGHT (LBS) 

S - AERODYNAMIC REFERENCE AREA < FT~2 ) 

ALGORITHM PARAMETERS: 

GL - GAIN FOR LIFT IN REFERENCE QBAR 

GHDOT - HDOT GAIN. EQ. GLIDE PHASE 
GHDOTX * HOOT GAIN. EXIT PHASE 
GO * QBAR GAIN. EQ. GLIDE PHASE 

GQX • QSAR GAIN. EXIT PHASE 
K 1 * DENSITY FILTER GAIN 

VTRIG * EXIT PHASE TRIGGER VELOCITY (FT/S^2) 

YBIAS * CONSTANT COMPONENT OF INCLINATION DEADBAND 

VSAT * DIVIDER FOR VARIABLE COMPONENT OF INCL. DEADBAND< FT/S ) 

SNKOFF * COSINE OF EO. GLIDE BANK ANGLE LIMIT 

SNKOFX* COSINE OF EXIT PHASE BANK ANGLE LIMIT 

GNC17B * CONSTRAINT FUNCTION GAIN 

QSTART * QBAR FOR GUIDANCE INITIALIZATION (LBS/FT~2) 

AMXER * INCLINATION DEADBAND 8IAS 
ZOFF » INCLINATION DEADBAND EPSILON 

RHONOT - ATMOSPHERIC OENSITY AT REFERENCE ALTITUDE (SL/FT~3> 
HNOT * REFERENCE ALTITUDE FOR EXPONENTIAL DENSITY (FT) 

HS • SCALE HEIGHT FOR EXPONENTIAL DENSITY (FT) 

RHOEX * OENSITY AT EXIT PLANE (SL/FT~3> 

GEXIT * GRAVITATIONAL ACCLERATION AT EXIT PLANE <FT/S^2> 

I GOX * NUMBER OF TIME STEPS TO WAIT BEFORE CALCULATING 
REFERENCE ALTITUDE RATE 
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DATA NP ASS 


DXYZIDI 3 > . 


VR. BANK. ALFBIA. 

XYZK3). XYZI D< 3 ) . 

. BANKCD. ALPCD , 

HEXIT. APOTG. DXYZK3) 

MU. RE. FPNM. RADCON. 

CL. CD. WT. S. 

GL. GHDOT. GHDOTX . GO. GOX. 

YB IAS . VSAT . BNKOFF . BNKOF X 
QSTART. AMXER . 20FF . 

RHONOT. HNOT . HS . RHOEX.GEXIT 
(RX XYZIUD. <RY. X YZI ( 2 ) ■> . <RZ. XY2I ( 3 ) ) 
( XDOT . X Y2I D t 1 > ) . ( YDOT . XYZID(2>>; (ZDOT. 

/ 1 / 


Kl. VTRIG. 
. GNC17B. 


I GOX 


XYZIDf 3 ) > 


C EXAMPLE VALUES FOR SELECTED INPUTS 
C 

C ■ PHYSICAL CONSTANTS 

C DATA RE / 20925741 . / . 

r • MU /I .4076469016/. 

C • FPNM /6076.U54/. 

C • RADCON /0. 0174532/ 

C 

0 VEHICLE PARAMETERS: 


C 

C 


DATA CL 

• CD 

• WT 

• S 


/. 29985/. 
/I .44164/. 
/ 224 1 7 . / . 

/ 1256 . / 


C ALGORITHM PARAMETERS: 
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CLBK 

C 


DATA 


DATA 


DATA 

DATA 


GL / -0 . 75/ . 

GHDOT /0 . 20/ . 

GHDOTX /0. 25/ • 

GQ / 5 . 0/ * 

GOX /0.0/. 

Kl /0 • 2/ . 

VTRIG / 23880 . / « 

Y B IAS /0 . 030/ * 

VSAT /26200./' 
BNKOFF /0. 997564/. 
BNKOF X /0. 965926/ 
GNC173 /3 . 6/ . 

QSTART /0.1/. 

AMXER /0 . 8 / . 

ZOFF /.01/. 

I GOX /2/ 

GEXIT / 3 1 . 0295/ . 
RHOEX /3.6164E-11/. 
HNOT /3 00000./. 
RHONOT /4.6235E-9/, 
KS / 1 7500 . / 


03/10/87 

IF( NPASS.GT. 1 )GO TO 2 


C 

C I N 1 T I L I 2AT IONS 
C 

G * MU / ( RE*RE ) 

M-WTZ/G 
I COUNT *0 
LATSV-0 
RLATSW*LATSW+1 
INHIBT-1 

CALL CROSS(OXYZIO.DXYZI.UYD) 
CALL UN I TV ( U YD . U YD ) 


NRR *0 . 


HDREF-0.0 
I GO* 1 
IEXIT-0 
KRHOD* 1 .0 
S 2ROL L » 1 .0 
WEDGM*0. 
WEDGA-0. 


SL P *0 . 
NPASS * 2 
RETURN 
C 

2 CONTINUE 
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RV « SORT ( RX*RX ♦ RY*RY * RZ*RZ> 

V » SORT! XDOT'XDOT * YDOT'YOOT + ZDOT*ZDOT > 
SGAM - DOT ( X YZ I . X YZ I D > / ( RV * V ) 

KDOT » V • SGAM 
VSQ -V/VSAT 
H • RV-RE 
GACT-MU/(-RV*RV> 


COMPUTE OSAR AS WOULD BE DERIVED * 
FROM PREDICTED CD 4 ACCELEROMETERS 


DRAGA ■ -XACC*COS( ALPHAR >-ZACC*SIN< ALPHAR > 
OPRED - DRAGA*M/< CD*S ) 


COMPUTE RHO BASED ON DRAG 


RHO - QPRED*2./VR**2. 

RH0S2 » RKONOT * EXP<< HNOT - H ) / HS ) 

CALCULATE RH0/RH062 AND FILTER TO CREATE 
DENSITY MULTIPLIER - KRHOD *•••••*• 

IF (V.LT.VTRIG) KRHOD»< 1-K1 > *KRHOD*K 1 *RH0/RH062 


CONTROL EQUATIONS 


......... CHECK QBAR FOR GUIDANCE INITIATION ••«••»»****•• 

I F < H . GT. 400000. )GO TO 100 

IF (QPRED .LT. QSTART ) GO TO 100 

QSTART-0. 

......... CALCULATE COS(PHI) FOR EQUILIBRIUM GLIDE •*»***• 

......... 2AN< ANGLE USED FOR VERT. LIFT MODULATION »****• 

COrHI » (WTZ/!CL*OPREO*S> > • (1.0 - V»*2 / ( GACT • RV ) ) 
COPHI1-COFHI 

....... Q 5A s reference equation for ecuil gl. »*«««***•* 

10 QREF*M*GACT/!CL*S*GL )•( 1 . -V**2/< GACT*RV ) ) 

* TEST FOR EXIT PHASE OR EQUIL. GLIDE •••*•*••**** 

IF (V.GT.VTRIG > GO TO 15 


• CALCULATE REFERENCE ALTITUDE RATE FOR EXIT PHASE **** 
.......... EVERY IGOX INTEGRATION STEPS *•*•***•««»•*»* 

ICOUNT » ICO'JNT ♦ 1 
I F ( ICOUNT. NE.IGOJGO TO 15 
CALL HDOTC< M. H .KRHOD .GAMEXD > 

• SET GAINS FOR GUIOANCE CONTROL EQUATION IN EXIT PHASE 

I COUNT *0 
IEXIT-1 
GO » GQX 
GHDOT - GHDOTX 
1GO-1GOX 

...... COMPUTE PLANAR WEDGE ANGLE *•«»•*•**•••*•••****« 


15 


CONTINUE 

SLPS-SLP 

WEDGAS-WED3A 

CALL CROSS! XYZID.XYZI .UYDS ) 
CALL UNITVf L'YDS .UYDS ) 
UYDOT*ABSt CCT( UYD.UYDS ) > 

I F ( UYDOT.GT. 1 . 1UYDOT-1 . 

WE D3A * ACOS ! UYDOT ) /RADCON 
SL?«W£DGA-WEOGAS 


CRIGSNAL PASH !S 

OF POOR QUALITY 
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*•**••* COMPUTE THE CONSTRAINT FUNCTION. Z 

FEQ*GNC17B«{ VSO-1 . ) 

IF(FEQ.LT.YBIAS/2. > FEQ- YB I AS/2 . 

Z-FEO-WEDGM 

IF(Z.GT.AMXER*VEDGM)Z«AMXER+VEDGM 
IF(Z.LT.YBIAS JZ-YBIAS 

IF< SLP-SLPS.GT.0. .OR. INHIBT.NE . I . OR . WEDGAS . GE .Z )G0 TO IB 
INHIBT-2 

IF{ FEO.GT. YBIAS/2. >G0 TO 17 
I F< VEDGAS.LT. UEDGM 1VEDGM-VEDGAS 
GO TO 18 

17 CONTINUE 

WEDGM-WEDGAS 

18 CONTINUE 

........ CHECK ROLL REVERSAL LOCK — * — — — 

I F ( INHIBT.EQ.2.AND.ABS(WEDGA-Z).GT.Z0FF HNHIBT-3 

......... SET RQLL REVERSAL command on 

IF(WEDGA.LE .Z.OR. I NH I BT . NE . 3 . OR . SLP .LE .0 . )GO TO Z0 
S2ROLL —S2ROL L 
NRR-NRR+1 . 

LATSW-0 

RLATSW-LATSV+1 

IF(NRR.EQ.2. )LATSW«2 

I F ( NRR . EQ . 2 . >RLATSV-LATSV+1 

INHI3T-1 


• — CALCULATE COMMANDED BANK ANGLE FOR EQUIL. GL. & EXIT *••••' 

.......... CALCULATE COMMANDED COSCPHI) — - — — * 

I F ( ( I EX IT .EQ. 1) .AND. (OPRED .LT. .5)) GHDOT- .tffftftfl 
ZB COPHI»COPHI-GHDOT-( HDOT-HDREF )/QPRED*GO»< QPRED-QREF > /OPRED 
COPHI2-COPHI 

I F < ( IEXIT.EQ.l ) .AN0.( OPRED.LT. .5 J.AND. ( IFIX.EQ. 1) JCOPHI-l .B 

**»* LIMIT BANK ANGLE TO 15-165 DEGREES BEFORE 2B0K FT. ***••*• 

I F ( IEXIT.EQ.l .AND . H . GT . 200000 . 1BNK0FF »BNKOFX 
I F ( COPH I .LT. -BNKOFF) COPHI — BNKOFF 
I F ( COP H I .GT. BNKOFF) COPHI- BNKOFF 

IF (COPHI . LT . -BNKOFF . AND . H . LE . 200000 . ) COPHI— BNKOFF 
IF (COPHI. GT. BNKOFF. AND. H. LE. 200000 . ) COPHI- BNKOFF 
IF (COPHI .LT. -BNKOFF. AND. H . GT . 3B00ET0 . ) COP H I --BNKOFF 
IF (COPHI. GT. BNKOFF .AND.H.GT. 200000 . ) COPHI- 8NK0FF 

PASS COMMANDED ANGLES BACK OUT TO ENVIRONMENT 

BANKCD- S 2 ROLL -AC OS ( COPHI >/RADCON 
ALPCD - ALFBIA 

100 CONTINUE 
RETURN 
END 
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c 

C 

C SUBROUTINE HDOTC 

c .................. 





SUBROUTINE HDOTC< M. H . KRHOO . GAMEXD ) 


c ........ Ex i T PHASE REFERENCE ALTITUDE ******** 

C **•***•*• RATE PREDICTOR/CORRECTOR ••••«»•• 

c .. ........ 

C 

C LOCAL VARIABLES 

C 

c 

C APOTG - TARGET APOGEE ALTITUDE <NM> 

C OIFF » DIFFERENCE BETWEEN PREDICTED APOGEE AND TARGET APOGEE (NM) 

C GAMEX - PREDICTED EXIT FLIGHT PATH ANGLE < RAO > 

C GEXIT - GRAVITATIONAL ACCELERATION AT THE EXIT PLANE <FPSS) 

C HAPO - PREDICTED APOGEE ALTITUDE (NM) 

C HOOEX » PREDICTED ALTITUDE ACCELERATION AT EXIT (FPSS) 

C HNOT » REFERENCE ALTITUDE USED IN EXPONENTIAL DENSITY MODEL (FT) 

C HS » SCALE HEIGHT USED IN EXPONENTIAL DENSITY MOOEL 

C QEX - PREDICTED DYNAMIC PRESSURE AT EXIT (PSF) 

C RE - EQUATORIAL RADIUS OF EARTH (FT) 

C REXIT » RAOIUS AT EXIT PLANE (FT) 

C RHOEXIT » ‘62 STANDARD DENSITY AT EXIT PLANE ( SL/FT3 ) 

C RHONOT ■ >62 STANDARD DENSITY AT REF. ALT. USED IN EXPON. MODEL ( SL/FT3 ) 

C VEXIT « PREDICTED EXIT INERTIAL VELOCITY (FPS) 

C VEX ITR - PREDICTED EXIT RELATIVE VELOCITY (FPS) 

C 


SAOD V: AOTVCM 

C 

C 

C IMPLICIT DOUBLE PRECISION (A-H. 0-2) 

REAL M.XRHOD 

C DOUBLE PRECISION M.MU . KRHOD . XI 

C 

C COMMON / HYPASCM / VR. BANK. ALFBIAS. DRAGA. FPACCI3). 

XY2I(3). XY2ID( 3 ) . 

BANKCD . ALPCO . 

HEXIT. APOTG. DXY2K3). DXY2ID(3), 

MU. RE. FPNM. RADCON. 

CL. CD. WT. S. 

GL. GHDOT. GHDOTX. GO. GOX. K1 . VTRIG. 
YBIAS. VSAT . BNKOFF . BNKOFX . GNC17B. 
QSTART . AMXER. 20FF. TSTEP . RRMAX . 
RHONOT. HNOT. HS. RHOEX. GEXIT 


V 

C • 

c * 

c 

c 

c * 

c 

c 

DELTA - 3 0. 

REXIT » HEXIT ♦ RE 
1-0 

C * S*CD/2./M 


C 

C ******* MULTIPLY DENSITY BY DENSITY SCALE FACTOR 
C 


RHONOT -RHONOT* KRHOD 


C 

c •«•«««• CALCULATE RHS OF PREDICTION EQUATION 
C 


IF ( HDREF . GT* 200. > GO TO 20 
HD RE F * 200 . 

20 RHS - C*RHONOT/HDREF*HS*( EXP( ( HNOT-HEX IT )/H$ )-EXP< ( 
♦ HNOT-H ) /HS ) ) 

1*1 + 1 
C 

C *********** CALCULATE VEXIT RELATIVE •*•**»•****•***•** 

C 


VEXITR - l./( l./VR-RHS) 

C 

C *********** BIAS VEXIT RELATIVE TO VEXIT INERTIAL 
C 


VEXIT - VEXITR + (V-VR) 


** ADD POTENTIAL/KINETIC ENERGY TRANSFER TO EXIT VELOCITIES ***** 

V2 ■ ( <GACT*H-GEXIT*HEXIT+(V»*2. )/2. >*2. )**0.5 

OELV - V2-V 

VEXIT - VEXIT OELV 

VEXITR » VEXITR ♦ DELV 
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HDOEX«VEXIT—2./REXIT-GEXIT-tCL-0.S*RHOEX*VEXITR-*2.*S>/M 

IF ( HDDEX.LT. 0. > GO TO 42 
IF < H.GT. 320000. ) GO TO 40 

DELHO « SQRT( HDREF**2 . *160000. *HD0EX ) - HOREF 
GO TO 41 

40 DELHD « SORTt H0REF«»2. *2. «H00EX»t 400000. -H > > - HDREF 
........... CALCULATE EXIT FLIGHT PATH ANGLE *••••***«**•***•**••*• 

41 IF{ (HOREF*DELHD ) .GT.VEXIT. OR. HDREF .LE.0. )G0 TO 56 
GAMEX « ASINttHDREF*OELHD ) /VEX IT) 

GO TO 45 

42 GAMEX a ASI Nt HOREF /VEX IT > 

........... CALCULATE APOGEE ALTITUOE ••••••••••***••*•***•*•****•• 

45 HAPO » (REXIT/(2.-REXIT*VEXIT**2./MU)*t l.*SQRT< I . REX IT*VEX IT* 
+COSt GAMEX ) )**2./MU*(2./REXIT-VEXIT**2. /MU ) > )-RE J/FPNM 


........... COMPARE PREDICTEO ANO OESIREO APOGEES ■**•***••*••** 

DIFF » HAPO - APOTG 
IF <ABS(OIFF>.LT.0.5) GO TO 56 

IF <SIGN( 1 .0*00. DIFF >.NE. SIGN! I. D*00.DIFFL) ) OELTA • DELTA/2. 
DIFFL ■ DIFF 

IF < HAPO. GT. APOTG) HDREF “HD REF -DELTA 
IF ( HAPO. LT. APOTG ) HOREF“HDREF*DELTA 
GAMEXO - GAMEX / RADCON 
GO TO 20 

S6 IFIHDREF.LT. 0. >HDREF«0. 

RETURN 

ENO 

SUBROUTINE UN ITV( A. B ) 

B • A / MAGNITUDE! A) 

IMPLICIT DOUBLE PRECISION (A-H. 0-2) 

DIMENSION A<3>. B(3> 

XMAG » SQRT< At 1 >*At 1 > ♦ A(2)*At2) * At3)»A(3) ) 

B< 1 ) - All) / XMAG 
B < 2 ) - At 2 ) / XMAG 
813) - At 3 > / XMAG 
RETURN 
ENO 

SUBROUTINE CROSS! A. B. C ) 

C ■ A X B 

IMPLICIT OOUBLE PRECISION 1A-H. 0-2) 

DIMENSION At 3 ) . BI3). CI3) 

Cll) » Al 2 ) • Bt 3 ) - A< 3 ) • B12) 

Cl 2 ) ■ At 3 ) * Btl) - At I ) * Bt 3 ) 

Ct 3 ) - All) • Bt 2 ) - At 2 ) * Btl) 

RETURN 

END 

OOUBLE PRECISION FUNCTION DOTt A. B ) 

RETURNS DOT PRODUCT OF A AND B 

IMPLICIT DOUBLE PRECISION tA-H. 0-2) 

DIMENSION At 3 > . 3(3) 

DOT - At!) * Btl) + At 2 ) • Bt2) ♦ At 3 ) • B<3) 

RETURN 

ENO 
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